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Abstract 

We detail numerical methods to compute the geometry of static vacuum black holes in 6 dimen- 
sional gravity compactified on a circle. We calculate properties of these Kaluza-Klein black holes for 
varying mass, while keeping the asymptotic compactification radius fixed. For increasing mass the 
horizon deforms to a prolate ellipsoid, and the geometry near the horizon and axis decompactifies. 
We are able to find solutions with horizon radii approximately equal to the asymptotic compactifica- 
tion radius. Having chosen 6-dimensions, we may compare these solutions to the non-uniform strings 
compactified on the same radius of circle found in previous numerical work. We find the black holes 
achieve larger masses and horizon volumes than the most non-uniform strings. This sheds doubt 
on whether these solution branches can merge via a topology changing solution. Further work is 
required to resolve whether there is a maximum mass for the black holes, or whether the mass can 
become arbitrarily large. 
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1 Introduction 



If one day we discover that there are extra dimensions in our universe, and these are well described 
by classical gravity, then they are likely (although not definitely pQ) to be compact along the lines of 
Kaluza-Klein theory Furthermore, if matter is confined to branes, then the radius of compactifi- 

cation could potentially be extremely large |U0)]. The simplest regular static vacuum solutions are then 
compactified uniform black strings [HlEj. Gregory and Laflamme (GL) showed that these are stable 
provided the horizon radius is large compared to the compactification scale |8|l9|ll0| . However, they also 
discovered a new family of non-uniform solutions emerging from the critical uniform string whose mass 
separates the more massive stable strings from the less massive unstable ones. These non-uniform solu- 
tions were constructed numerically, firstly by Gubser as a perturbation expansion in a non-uniformity 
parameter A about the A = critical uniform solution jllj . and then non-perturbatively in |12| using 
elliptic methods. The third class of solutions expected to exist are black holes that do not wrap the 
circle direction. In 4 dimensions such solutions were found by Myers analytically ( see a ls° gener- 
alisations jl4U15j with modified asymptotics) but in more than 4 dimensions little is known |16U17j . 
essentially as the rotation group then has curvature. Using the elliptic numerical methods of |18lll2j 5 
dimensional localised black holes have recently been constructed on a Randall-Sundrum brane |19U2U| . 
This is a related numerical problem of considerable interest as recent conjectures claim large localised 
static black holes may not exist |21[ I22|. 

Kol proposed an elegant relation between these three types of Kaluza-Klein solution [2Sj, the non- 
uniform strings linking the uniform branch to the black hole branch. The string to black hole transition, 
also explored in |17 [l24H25| l2fi | l27 | l28| . is then conjectured to be continuous, and have a Lorentzian cone 
geometry where the horizon degenerates, and this agrees very well with numerical tests on the non- 
uniform string branch |29U3U| . If this picture is correct, it predicts that the black holes, like the 
non-uniform strings, have a maximum mass. While for fixed compactification radius the 3 classes of 
solution overlap at intermediate mass scales, 1 at large masses the uniform strings would be the unique 
non-singular solutions in Kaluza-Klein theory. If incorrect, it may be possible to have arbitrarily high 
mass black holes, if either the geometry becomes increasingly 'squashed' or decompactifies on the 
symmetry axis (as in pure Kaluza-Klein theory there is no radius stabilisation). Whether this could 
persist in a radius stabilised theory would then be an important phenomenological question. 

Clearly Kaluza-Klein theory is a simplification of realistic compactifications. Kol's picture presum- 
ably remains unchanged adding warping |33H 34 |l331 l3fi|. or charging black holes under matter fields 
localised to branes jSZ]. For additional bulk matter, such as is necessary for stabilisation, the situation 
may be more interesting, but we expect it will inherit many features of the pure Kaluza-Klein case. 

The Gregory-Laflamme instability, underlying the dynamics of these compactified horizons has been 
linked to thermodynamic stability [3*51I391I4U1I4HI42U431I36U441I45| , and an analogous classical instability 
has recently been conjectured for the rotating Myers- Perry solution |46U47lH5j . which is thought to 
be unstable for large angular momenta. The end-state of the classical Gregory-Laflamme instability is 
still a mystery, although there has been interesting analytic and numerical work on this subject |49|I50| . 
We note that stable black strings evaporating via Hawking radiation will eventually succumb to this 
classical instability, and understanding the dynamics, and in particular whether cosmic censorship is 

x We must consider the full asymptotic charges to distinguish the solution |26U27| . although whether there is a unique 
solution with these charges is an interesting open question |28|. as uniqueness constraints apparently weaken in more than 
4 dimensions 
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violated, is important in order to understand evaporation of cosmological/astrophysical black holes 
below the compactification mass scale. 

The objective of this paper is to numerically construct and study the non- wrapping black hole 
branch of solutions. We perform this analysis in 6-dimensions so that we may compare with the 
previous non-uniform string numerical results of ^21 > which for technical reasons, were performed 
in this number of dimensions. We begin with a brief discussion of the numerical method, which 
involves phrasing a subset of the Einstein equations in a way compatible with numerical relaxation, and 
most importantly, showing how the remaining 'constraint' equations can be satisfied by appropriately 
choosing the boundary conditions. Since the method has now been used several times |18|I12| IT9] we 
refrain from a detailed exposition, and instead highlight the various subtleties related to this Kaluza- 
Klein black hole problem. We then go on to discuss the numerical results. We demonstrate that as 
expected these solutions exist, at least within the scope of our numerical approximation. 

We compute geometric embeddings of the spatial horizon and symmetry axis into Euclidean space, 
and show that the geometry near the axis decompactifies with increasing mass, and the event horizon 
deforms to a prolate ellipsoid. With the current implementation we are unable to ascertain whether 
this decompactification terminates with a maximum mass black hole that just 'fits' into the compact 
direction, or whether the decompactification continues indefinitely so that arbitrarily high masses can 
be found. 

The maximum size black holes we are able to construct have horizon radii approximately equal to the 
asymptotic compactification radius. We compare these with the most non-uniform strings constructed 
in ^21 finding the mass and horizon volume of these modest sized black holes already becomes larger 
than that of the maximally non-uniform strings, and the axis decompactifies to a greater extent. The 
size of black hole we may construct is limited by numerical factors, and it seems clear that still larger 
black holes exist, with the above trends continuing for these. The implication is that it appears unlikely 
that the non-uniform string branch (connected to the critical uniform string) and this black hole branch 
are connected via a topology changing solution. 

Various technical details and numerical checks are reserved for the three Appendices. We pay 
particular attention to ensuring and checking that the constraint equations are indeed satisfied for the 
solutions. 

The reader is also referred to independent work by Kol, Piran and Sorkin who we understand have 
recently performed related calculations in 5 dimensions |51| . 

2 Method 

In order to solve the black hole geometry we are required to solve the Einstein equations with elliptic 
boundary data; we wish to have a regular horizon geometry, for the solution to be periodic, and also 
asymptotically to tend to flat space product with a circle. We employ the methods first developed in ^B] 
and used to construct non-uniform string solutions ^2] , and later localised black holes on branes [T31I2U] . 
2 In this section we outline the method and boundary conditions appropriate for the problem. Due to 
its necessarily technical nature, some readers may wish to skip to the following 'Results' section. For 
a more general discussion of the method, the reader is referred to 12 . Technical numerical details are 

2 See also |52| for a method of solving static axisymmetric non-vacuum black holes in 4-dimensions which shares some 
features with our method. 
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also provided for the interested reader in Appendix^ and important numerical checks are reported in 
Appendix El to demonstrate the method performs correctly. 

Constructing the non-uniform strings is a very clean situation in which to apply these elliptic 
numerical methods. However, the black hole problem at hand is substantially more difficult, primarily 
for two reasons; 

Firstly weakly non-uniform strings can be described as a perturbative deformation of the critical 
uniform string. The relaxation methods employed here require a good initial guess, or typically no 
solution will be found. Thus for the strings the non-uniformity can be turned on 'gently'. In analogy, 
a very small black hole will appear as a 6-d Schwarzschild solution near its horizon, but it obviously 
must have very different asymptotics due to the compactification. Thus even for a small black hole, we 
do not have an exact solution to 'gently' start building larger black holes from. We tackle this issue by 
building in 6-d Schwarzschild behaviour at the horizon that decays quickly away far from the horizon, 
and then we solve for the correction to this, which should be small for low mass black holes. 

Secondly, the axis of symmetry is exposed in the problem. The coordinate singularity at the axis 
generically gives rise to problems numerically and there are various ways around this in conventional 
evolution problems (see for example |53j). As discussed in |18j the elliptic method we use is very 
sensitive to this coordinate singularity, which may destroy the ability of the algorithm to relax to the 
solution. Furthermore the coordinate system we require to phrase a subset of the Einstein equations 
in an elliptic manner introduces even worse coordinate problems on the symmetry axis than one would 
normally expect. We have found no elegant method to tackle this problem, but do have a functional 
approach, originally used in jlSj and discussed here in Appendix Improving or evading this problem 
appears to be crucial for increasing the capability of this method. 

So whilst the problem is a rather delicate one, we are still able to make progress. As in previous 
applications of the method, we write the static axisymmetric metric in a diagonal form, retaining a 
conformal invariance in the radial and tangential coordinate r, z as; 

ds 2 = g MN dx M dx N = -e 2a dt 2 + e 2 ^ {dr 2 + dz 2 ) + r 2 e 2(3+ hdn 2 3) (1) 

with a, /3, 7 being functions of r, z. The particular linear combination of 0, 7 taken above is simply for 
technical convenience later. We take the z coordinate to be compact with period L, and later will require 
the metric functions a, (3, 7 to vanish at large r, and hence the physical radius of compactification will 
be L. We choose units such that the 6-dimensional Newton constant, Gjv(6) = 1- Since we may perform 
a global scaling on any solution of the vacuum Einstein equations, for future convenience we chose to 
set L = 7r in these units. 

One nice property of this form of the metric is that one can choose the position of the boundaries 
in the (r, z) plane to be at any location, due to the residual conformal coordinate transformations. A 
second important feature of this coordinate system is that 3 of the 5 Einstein equations, G t t , (G p p + 
G x x ),G e g have elliptic second derivatives, being just the (r, z) Laplace operators, and thus we term 
these the 'elliptic' equations. As we are so far unable to write a positive definite functional of the 
metric components which can be minimised to give these equations, it is not at all clear the problem 
is truly elliptic. However one can still use relaxation methods to solve them, specifying elliptic data 
on the boundaries of the problem. The most important feature of this coordinate system is that using 
the contracted Bianchi identities, the 2 remaining 'constraint' equations weighted by det gu n, 

$ = det g MN G r z # = l - det g MN {G r r - G\) (2) 
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obey Cauchy-Riemann (CR) relations, 



«9 r $ = d z i£ 



(3) 



if the elliptic equations are satisfied. Simply relaxing the 3 elliptic equations for the 3 metric functions 
will generically yield a solution, but this is only consistent with the full set of Einstein equations if 
the boundary data is such that the CR relations imply both constraints are satisfied. For example, for 
the non-uniform strings we may impose the constraints by updating the elliptic boundary data such 
that the G r z constraint is satisfied on all boundaries, and the remaining constraint (G r r — G z z ) is just 
imposed at one point, solving the CR problem )12|. In our example here, we in fact find it convenient 
to impose both constraints, but on different boundaries, in such a way that still provides sufficient 
conditions to satisfy the CR problem, but does not over determine the elliptic equations. 

2.1 Boundaries and coordinates 

We still have residual coordinate freedom, and we use this to tailor the coordinate system to our 
problem. Instead of relaxing a, /3, 7 directly, we wish to perturb about a 6-d Schwarzschild solution 
near the horizon, 



where the 'background functions' {A, B,C}b g suitably express this Schwarzschild geometry at the 
horizon, decay away radially, and are compatible with the compact boundary conditions - so we could 
obviously not just take the Schwarzschild metric itself. Notice that C = corresponds to conformally 
flat spatial sections. 

Let us now consider a form for the background functions, and boundary locations in the (r, z) 
coordinates. The 6-d Schwarzschild metric can easily be written in an appropriate conformal form as, 



with fj? = r 2 + z 2 . To be suitable for the background functions, this must be modified away from the 
horizon to ensure compatibility with the periodic boundary conditions in z. However, assuming after 
this modification the horizon remains at constant [i = fio, it will form a circle in the (r, z) coordinates. 
Thus we would wish to take boundaries of the form in figure ^ in order to represent the coordinate 
axis and black hole horizon, and periodic boundaries. Then we would reasonably expect the functions 
A, B, C to remain finite everywhere, and be small for a small black hole, allowing us to use the initial 
data A = B = C = for the relaxation. 

Numerically it is always convenient to have a rectangular grid. Whilst we are in principle free 
to use the residual conformal transformation to fix the boundaries to be wherever we wish, clearly 
to obtain a rectangular domain such a transformation must be singular, as one right-angle in the 
figure ^ should be 'flattened' out. However, if we find such a coordinate transformation analytically, 
we may separate out any singular behaviour from A,B,C, leaving their behaviour perfectly regular. 
Any conformal coordinate transformation is generated by a solution to the 2-d Laplace equation, and 
choosing a solution p(r, z) to be that representing a point source in the compact 2-d space (ie. on a 



a = A bg + A (3 = B bg + B 1 = C bg + C 



(4) 
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Figure 1: Schematic illustration of the boundaries we would intuitively take in r,z coordinates. At 
least for small black holes, taking the horizon boundary to be circular would ensure the metric functions 
A,B,C would be small and finite, as the geometry near the horizon would only be a weak distortion 
of 6-dimensional Schwarzschild, which we build into the metric ansatz via A^ g , B^g, C\, g . 



cylinder) we may define new coordinates (p, x), 



p(r, z) 
X(r, z) 



- log 

2 & 



tan 



— (cosh 2 r — cos 2 z) 



tan z 



(6) 



tanhr 



where x is determined from p by CR relations. These essentially 'flatten' out the horizon, and now 
X is the compact coordinate which conveniently takes the range of an angular coordinate (0, 7r/2) for 
half a period of the solution, z = (0,L/2) = (0,7r/2). We illustrate the isosurfaces of p and x m the 
(r, z) plane in figure |2j Contours of constant p, for p < 0, generate very similar curves in the (r, z) 
plane to that of the horizon in figure ^ Note also that x = 7r /2 gives us both the axis of symmetry (for 
p < 0) and the periodic boundary z = ir/2 (for p > 0), and p = 0, x = 7T /2 is the singular point in the 
conformal transformation. Thus if we use these coordinates, a rectangle with x = to ir/2 and p = po 
to oo will, for po < 0, give us similar looking boundaries to figure ^ 

Whilst the coordinate transform is rather singular, since we have an analytic expression for it, we 
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Figure 2: Illustration of isosurfaces of p, x as functions of r, z. The schematic boundaries of figure 
Q can be mapped to constant p and x values, making for convenient numerical implementation. For 
p < 0, x behaves as an angular coordinate, and p as a radial one, whereas for large positive p, we find 
p,X behave as r,z respectively. We will take the horizon boundary at constant p = po, finding later 
that specifying po determines the physical size of the black hole solution. 



may remove the singular Jacobian, J, and now write the metric, 



ds< 



-e 2a dt 2 



J(p,x) 



r(p,X) 



+ 



^-^J{p, X ) (dp 2 +d X 2 )+r(p,x) 2 e 2 ^hdnl 

P 2p 



(3) 



(7) 



y/1 + e 4 P + 2e 2 P cos 2x 



1 



■ cosh 



-i 



=2p 



+ ^l + e^P + 2e 2 P cos 2 X 



where now a, f3, 7 are exactly the same functions as previously in (^Q), except now in terms of p, x- 

Now consider the 6-d Schwarzschild metric ®. For small po, near the horizon we find e 2p ~ p 2 = 
r 2 + z 2 , with p behaving as a polar radial coordinate and x as t ne angular one. Since we wish to choose 
background functions to reproduce the Schwarzschild metric near the horizon for small black holes (ie. 
small po), and yet to implement the compactness requirement in z, a simple choice is just to substitute 



0P-P0 



p/po into (JSJ) giving, 



e B bg 



1 



=3(p -p) 



1 _|_ e 3(po-p) 



(8) 



1 + e 



3(po-p) 3 



Cbg = 

The independence of these background functions on x ensures they satisfy the periodicity requirements. 
The constant po tells us the coordinate position of the horizon. We will shortly show that keeping the 
asymptotic compactification radius fixed, it is the parameter p$, and thus the contour p = po we take 
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to be the horizon that we will vary to change the mass of the black hole. From the earlier figure |2] 
we see that po must be negative to give a spherical topology for the horizon. The small horizon limit 
is now po ~~ > —oo. For po < but closer to zero, the ansatz above gives a deformed geometry from 
Schwarzschild too. 3 Since we include the vanishing of the lapse in A^ g at p = po we expect the metric 
functions, and in particular A, the one associated with deformations of the lapse, to remain finite 
there. Far away from the horizon, ie. for p — » oo or alternatively r — ► oo, the background functions 
decay exponentially in p or r. We stress that for finite negative po, A, B, C = is not a solution to 
the Einstein equations, but for very negative p, the small black hole limit, A,B,C will at least be 
small everywhere as the horizon tends to Schwarzschild, and by the time the metric functions 'see' the 
compactification they will be vanishingly close to zero anyway. After the change in coordinates, we 
obtain analogous elliptic equations, G t t , (G p p + G\), G e and also CR relations for the new constraint 
functions <£, fy, 

$ = det qmnG^ * = X - det g MN (G% - G\) (9) 

so that, 

d p $ = d x V d x $ = -d p ^ (10) 

2.2 Boundary conditions from the elliptic equations 

The 3 elliptic equations we are solving require various boundary conditions due to the regular singular 
or periodic behaviour at these boundaries. To satisfy the constraint equations we must impose more 
than just these conditions. However, let us start by considering the basic boundary conditions from 
the elliptic equations. 

Asymptotically we want the geometry to be a product of 5-dimensional flat space with a circle, and 
thus we take A,B,C — > 0. Since we earlier fixed the range of %, this also fixes the compactification 
radius to be L = it. Of course we may simply globally scale these vacuum solutions to obtain any 
desired asymptotic radius. We find the asymptotic form required by the 3 elliptic equations, 

A~*+o(\) B ~*+± + o(\) C ~ h + *L + (\) (11) 

p 2 Vp / p p \p J p p \p J 

25 

with 3a 2 + 96 2 -4c 2 H b\ = 

6 

expanding (without linearising in the metric components) in inverse powers of p ~ r. Fourier modes 
with x dependence decay exponentially (since x — z i f° r large p), as does the contribution from 
the background functions {A, B,C}b g - Thus on our asymptotic boundary we have mixed Neumann- 
Dirichlet boundary conditions for the 3 metric functions. 4 

The symmetry axis x = 71 /2 with p < requires that A,B be even in (x ~ 7r /2) (or alternatively 
r). However we also find the requirement that C = as there is a regular singular behaviour due to 

3 It is interesting to compare this coordinate system with that proposed by Harmark and Obers in |17| where the radial 
coordinate used followed a 6-d equipotential. Here, our coordinate follows a 2-d equipotential. Presumably using the 
higher dimensional potential is a sensible procedure, and could give an improved ansatz to perturb about, which better 
models the horizon geometry for the larger deformed black holes. It remains an interesting problem where one can use 
the ansatz of Harmark and Obers to do numerics with, particularly as the ansatz and having the horizon at constant 
potential location was proven to be consistent for non-uniform strings in 129) and recently for black holes in |28j 

4 In practice we impose these conditions at a finite, but large p — p ma x and check in Appendix^] that the results are 
independent of pmax- 
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the form of the coordinate system. We might be confused that there is no Neumann condition on C 
but we see later this emerges from the constraints. 

We require that the metric functions A,B,C be finite at the horizon p = pq. With the form of 
background functions {^4, B, C}b g we then find that the elliptic equations are consistent provided that, 

Ap\p=po = (12) 

r _i_ _r l — i _ -£ l 

)P 3 \p = pq — 1 \p=P0 

with r(p, x) as given earlier in equation (J7J). 

Finally, the elliptic equations at the remaining periodic boundaries at x = an d X = 7r /2 with 
p > simply imply Neumann conditions on the metric functions A, B, C. 

2.3 Boundary conditions from the constraints 

The constraint equations also impose conditions on the metric. We have seen that assuming we can 
satisfy the elliptic equations for given boundary data, the two constraints obey CR relations. Thus we 
do not need to enforce both constraints on all the boundaries, which naively would over-determine the 
elliptic equations. Firstly we consider the extra conditions the constraints impose, and then discuss how 
best to implement them to ensure a consistent solution of the CR problem, without over determining 
the elliptic data. 

On the symmetry axis and periodic boundaries we have already specified 3 conditions, one for each 
metric function, and consequently treating this as a boundary value problem, we do not wish to impose 
any more. On the x = periodic boundary G p x vanishes by symmetry, and consequently so does the 
corresponding weighted constraint <£, and the remaining constraint {G p p — G\) is guaranteed to be 
even. A similar situation occurs on the x = 7r /2 boundary for p > which represents the other periodic 
boundary. Indeed these periodic boundaries are fictitious in the sense that we can consider the problem 
on the unwrapped covering space where these boundaries are 'removed', and thus we should not need 
to impose any constraints here. 

For x = 7r /2, but now with p < we have the symmetry axis. Again this is in principle a 'fictitious' 
boundary, but we must impose C = here for the elliptic equations, and it is hard to see how this 
would lift to a covering space with 'no' boundary. Thus we examine the situation at the axis in more 
detail. Using the boundary conditions from the elliptic equations, and assuming the metric components 
are regular, G p x vanishes there, but {G p p — G\) does not unless the normal gradient of C vanishes. 
However, since the measure det gMN vanishes near this axis, both <3? and are zero, and so from the 
point of view of solving the CR constraint equations we need do no more here. As discussed in the 
original implementation of this method |18j . this resolves the paradox that C = and C has a Neumann 
boundary condition, despite us solving C using a boundary value formulation. Whilst we only impose 
C = to be compatible with the elliptic equation behaviour, if we consistently provide data on the 
other boundaries to ensure the weighted constraints vanish everywhere, then this Neumann condition 
follows automatically, and does not need to be explicitly imposed at the symmetry axis. 

Now let us consider the remaining boundaries, the horizon and asymptotic boundary. At the horizon 
we only have 2 conditions from the elliptic equations and require another to specify elliptic data. Firstly 
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at the horizon we find that the G p x constraint implies the horizon temperature is a constant, 

{A-B + C)J p=P0 = l -^f\ p=P0 (13) 
and the remaining constraint (G p p — G\) requires, 

(A-B + C) jP \ p=P0 =^\ p=P0 -l (14) 

Lastly, asymptotically at large p (or equivalently r) due to the exponential decay of x (or equivalently 
z) dependence the G p x constraint goes exponentially to zero, and so the weighted constraint $ is 
guaranteed to go to zero, even though detgMN ~ P 3 - Due to this power law growth of the measure 
det gMN, it is less obvious the (G p p — G\) constraint weighted as $ goes to zero, as (G p p — G\) 
(unlike G p x ) depends on the homogeneous components of the metric, which only decay as a power law. 
However, the behaviour of the homogeneous component implied by the elliptic equations in (jllj) does 
also ensure that ^ = asymptotically. 

Now we must decide how to specify data for the elliptic equations (ie. one condition for each metric 
function on each portion of boundary), but also satisfy the constraint problem. With the conditions 
already required by the elliptic equations, we have sufficient data and $ = on all boundaries except 
the horizon. The second weighted constraint \& = is satisfied only asymptotically. At the horizon 
neither constraint is satisfied, and we require one more condition for a linear combination of A,B,C 
to make up the elliptic data as we so far only have (|12j) . 

We could impose the constant horizon temperature condition and thus set $ to zero on the horizon. 
As we have seen $ is zero on all the other boundaries, and from the CR relations it obeys a Laplace 
equation, so this would uniquely set it to zero everywhere. ^ would then be zero following from the 
CR relations and the fact it vanishes asymptotically. 

However it is numerically more stable to impose the (G p p — G\) constraint, and hence ^ = 
at the horizon instead. The {G p p — G\) constraint is a typical elliptic boundary condition, whereas 
imposing the constant horizon temperature G p x constraint involves 'less local' tangential derivatives 
on the horizon. Now we have sufficient data to impose both constraints globally via the CR problem. 
Since at the horizon \E r = 0, $ has a Neumann condition there, and is zero on all other boundaries, and 
hence will be zero everywhere as it obeys a Laplace equation. The CR relations then imply that is 
constant, and must be zero as it was imposed to be zero on the horizon and is also true asymptotically. 

2.4 How p specifies the size of the black hole 

Fixing A,B,C — > asymptotically, and thus the asymptotic compactification radius to be L = tt, 
there must be one constant entering the boundary data that specifies the size, or mass of the black 
hole. Intuitively one would imagine this to be po as this certainly enters into the boundary data at the 
horizon, in equations (|12|) and (|14j) . 

To confirm this, we must demonstrate that po is a physical quantity, and not simply a coordinate 
artifact. Thus we must show there is no residual coordinate transformation (ie. conformal transforma- 
tion on p, x) that preserves the rectangular boundaries, and conditions on these boundaries, and the 
asymptotic radius of compactification, but changes the effective po in the new coordinates. If this were 
the case, po would not correspond to a physical parameter, and therefore could not specify the size of 
the black hole, which certainly is a physical parameter. 
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Let us suppose we have a black hole solution with a particular p . Now let us construct the most 
general coordinate transformation that simply preserves the rectangular boundaries. We construct 
new coordinates p — ► p = p(p,x), X ~~ * X = X.{PiX)i an d then must solve a CR problem to build 
the conformal coordinate transformation. Let us do this by specifying data for the Laplace equation 
that determines x.- F° r the boundaries to remain rectangular, we must specify constant x = Xo on 
the x = boundary and constant x = X-k/2 on the x = 7r /2 boundary. At the horizon, we must 
specify a Neumann condition on x an d then the CR equations guarantee p is constant there. Regular 
asymptotics then give the unique solution, 

X = X0 + (Xvr/2 - Xo) —i^ ( 15 ) 
' 7T/Z 

This solution now completely determines p up to a further constant of integration. 

Now that we have the general transform preserving the rectangular boundaries, let us further 
restrict it by making it preserve our boundary conditions. Firstly, since xo does not enter into any 
equation or boundary condition, we may freely set this to zero. Secondly, since we have selected 
A,B,C —* asymptotically, we require x = n/2 on the x = 7r /2 boundary if we are not to change the 
compactification radius L = it. Thus now, x = X- This implies from the CR relations, 

p = 5p + p (16) 

where 5po is a constant of integration. 

However there is a subtle point. Whilst po only enters the boundary conditions explicitly at the 
horizon, the boundaries conditions change on the x = 7r /2 axis for p > (where C has a Neumann 
condition imposed) to p < (where C = is imposed). Since we have fixed this transition of boundary 
conditions to occur at p = 0, in the new coordinates this will occur at po = 5po- Hence the transformed 
solution will not satisfy 5 the boundary conditions in the new coordinates where this transition now 
occurs at p = 0. Thus, due to our fixing the transition from symmetry axis to periodic boundary at 
p = 0, there are no residual coordinate transformations mapping solutions with horizon position po to 
a transformed solution, solving the boundary conditions but with a different horizon position p~o in the 
new coordinates. Therefore po does indeed specify the physical size of the black hole. 6 

5 Note that if the physical solution had C — and C, x = on all of the x = 7r /2 boundary then po would not specify 
the solution. As discussed the constraints ensure that whilst we only impose C = 0, C also satisfies a Neumann condition 
on the symmetry axis, p < and \ = t/2. However, for p > on the periodic boundary we only impose a Neumann 
condition and there is no reason why C would vanish there too, and indeed in the solutions we find C does not vanish 
there, although it is very small. 

6 It is interesting to note that the horizon must have spherical topology for this argument to work. If we were considering 
a string horizon, and thus chose po > 0, then the boundary conditions would simply be Neumann all along the \ ~ 
boundary, and the new transformed solution in the p, \ coordinates would successfully satisfy all our boundary conditions, 
but with a different value of po po) on the horizon. It is for this rather subtle reason that the boundary conditions in 
the non-uniform string case of II 21 must be imposed differently, using the G p x constraint equation on the horizon (rather 
than (G p p — G x x ) as we use here) which, being a tangential condition, rather than a condition on the normal derivatives, 
introduces a new integration constant that parameterises the mass, or equivalently A for the string solution. The fact 
that this is a tangential condition appears to make the algorithm 'less local' and to require considerable under-relaxation, 
whereas here we do not need this. However, here while we need not damp the relaxation, the exposed p < symmetry 
axis does lead to the coordinate singularity induced stability problems discussed above. 
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2.5 Thermodynamic quantities 

We will compute the temperature T and entropy S of the black hole solutions, and these are given as, 

1 Ze A - B+c , 



T 



27T 28/3^J ' P=P0 

vr/2 



(17) 

P=PO 



The mass may be computed by two independent methods. Firstly we may determine the mass from 
the asymptotics of the metric, as, 

ir 2 L 

M = — (-3oa - b 2 + c 2 ) , (18) 

using the expansions pip |54j . Secondly the mass may be determined by integration from the First Law, 
using dM/dpo = T(p )dS /dpo (which applies for fixed asymptotic compactification radius - see below) 
along the branch of solutions, and taking M = for po — ► — oo to define the integration constant. Later 
(figure [SJ we will see very good agreement between these two values. This is a good indication that 
the elliptic equations are well satisfied globally and the boundary conditions are imposed correctly. An 
important point is that the First Law does not test whether the constraint equations G r z , (G r r — G z z ) 
are satisfied, and we elaborate on this point in Appendix [UJ Since we completely relax the elliptic 
equations, it is then not terribly surprising that the First Law holds very well. What is absolutely 
essential is that the constraints are also checked, to ensure they are well satisfied. They indeed are, 
and these tests are outlined in Appendix IbI 

Since the black hole geometries are compactified, the mass is not the only asymptotic charge. This 
feature, shared more generally by branes was originally discussed in |55j . and more recently in relation 
to the black hole/black string problem in |26 |I27|I28| . It is easy to see there must be another charge. As 
the geometry becomes homogeneous at large distances from the symmetry axis it can be dimensionally 
reduced to Einstein-dilaton-Maxwell theory. For our regular static solutions the Maxwell vector can 
always be gauged away, but we are left with gravity and the dilaton scalar, indicating we should consider 
both an asymptotic mass, but also a scalar charge. 

From a purely 6-dimensional point of view, the second charge can be thought of as a binding energy 
per unit mass, n, resulting in a modified First Law, 

dM = TdS + nM^- (19) 

with the new term representing work done when varying the asymptotic size of the extra dimension. 
In our solutions, fixing L we reproduce the usual form of the First Law for black holes. However, we 
may determine n from the asymptotics of the metric, or from the Smarr relation, 

X ft 

TS = — M (20) 

again discussed in [HE], an d given explicitly for the problem at hand in 26,27 , where n was calculated 
from the asymptotics to be, 

a 2 + 36 2 -3c 2 . . 

n = "5 ZT ( 21 ) 

3 a 2 + 6 2 - c 2 
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As emphasised in |27j we may use the Smarr formula as a check of our numerics. This is very similar 
to the First Law check discussed above, which only involves the mass, whereas Smarr's law involves 
both M and n. It is important to note that for the same reasons the First Law does not probe the 
constraint equations, the Smarr formula also does not. Thus again it is no replacement for the checking 
of constraint equation violations we perform in Appendix IE! 

2.6 Behaviour of the method 

Following the above method, we may construct a unique numerical solution for each value of po. The 
elliptic equations can be solved very stably, once the coordinate induced instability at the axis is dealt 
with (see Appendix 0). Since we have designed the coordinates and method to have A, B,C — ► in 
the small black hole limit po — > — oo, the method behaves well in this limit. As the black hole becomes 
larger, so po is finite and negative, A,B,C deviate away from zero, although remain regular as we 
expect. Topology of the p coordinate dictates po = is the largest black hole that could exist. Using 
reasonable resolutions (up to ~ 140*420), we were able to find solutions with a maximum po — —0.18, 
yielding a black hole horizon with typical radius comparable to the asymptotic compactification radius. 
As an example we show the metric functions A, B, C for po = —0.28 in figure|3J Note that the magnitude 
of C is much less than that of A, B, and this is increasingly true the smaller the black hole. Hence the 
spatial sections are approximately conformally flat. 

The larger the black hole is, the larger the gradients in the metric functions, and for a fixed 
resolution the method no longer converges past a certain black hole size. Going to a higher resolution 
we find the problem is removed, and the size can be further increased, but obviously the problem then 
re-occurs at a new larger size. The key area where we lack resolution is near the symmetry axis. For 
the large black holes, with po closer to zero, there become fewer and fewer points there. With the 
maximum size we could find, po = —0.18, so the coordinate distance of the symmetry axis is 0.18, 
compared to the coordinate distance along the horizon which is ir/2 ~ 1.6, and hence with our simple 
discretization scheme (see Appendix \K\ for details) the axis is allocated far fewer points. The closer po 
gets to zero, the more acute this problem. Thus in our simple numerical implementation, we are limited 
by resolution, and hence computation time. We present results in this paper using modest resources 
and simple relaxation algorithms. It is likely that with improvements in both areas one can achieve 
far improved data. For example, adaptive grid methods may circumvent the lack of resolution near 
the symmetry axis, but it is a serious challenge to implement these and maintain a stable relaxation 
agorithm. However, already with our simple implementation it is possible to derive interesting physical 
results as we shall see. 

3 Results 

The questions we wish to address are whether there is an upper mass limit for these solutions, and 
whether the geometry is compatible with continuation to the non- uniform string branch. Whilst po = 
is the largest black hole due to the topology of the (p, x) coordinate system, obviously the metric 
functions A,B,C may diverge in this limit, and consequently so might the horizon volume and mass, 
and thus a priori we have no reason to assume such an upper mass limit will exist. On physical grounds 
one could argue that a black hole would not be able to 'fit' into the compact extra dimension, but we 
stress that in the unstabilised pure Kaluza-Klein theory it is only the asymptotic radius we have fixed 
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Figure 3: Plots of the metric functions A, B, C for a black hole solution with p$ = —0.28. This is 
quite near to the maximal size (po = —0.18) we were able to construct before we become limited by 
gradients and lattice resolution near the symmetry axis at x = 7r/2 for p < 0. Already this black hole 
solution has equal horizon volume and mass to the most non-uniform strings compactified on the same 
asymptotic radius. Whilst the lattice is large, being 140 * 420 in x, p, since po is very close to zero 
the number of points along the symmetry axis is only around ~ 20. This is still enough to see good 
behaviour in the metric functions. Note that B, C are much less than A in magnitude. The maximum 
p for the lattice is ~ 5, and not all the domain is shown as the functions simply go smoothly to zero at 
large p. 
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for the branch of solutions, and there is nothing to prevent the geometry along the axis and horizon 
from decompactifying as the black hole becomes larger. 



3.1 Horizon geometry 

In order to approach these questions we embed the spatial horizon geometry into 5-dimensional Eu- 
clidean space, 

ds 2 = dX 2 + dY 2 + Y 2 dQ 2 {3) (22) 
and matching our geometry at p = po implies that, 

Y(x) = re^h\ p=po (23) 



X( X ) = f X d X \ \j(x')e 2 ^^-p 2 

Jo V L '* . p=po 

and we interpolate the numerical data to perform the integral. Clearly for small black holes we expect 
a spherical horizon, and for larger black hole we expect deformation. We find excellent agreement with 
the horizon being a prolate ellipsoid for all our solutions up to the maximum size available po = —0.18. 
In figure 0] we plot a moderate and a large black hole to demonstrate the accuracy of the ellipsoid fit. 
We plot the positions of the actual lattice points in the embedding coordinates X, Y for our highest 
resolution, and against these we plot the fitting ellipse, and note that all the points fall consistently on 
the fit curve. Thus from now on it is easier for us to characterise the geometry using the major (polar) 
and minor (equatorial) axis radii, which we term R po i ar , R eq . Then using this elliptical fit we plot the 
ellipse radii and ellipticity, 

e = (24) 

against po for all our solutions in figureEl We see that for the largest black holes the ellipticity decreases 
to only ~ 0.87, even though the ellipsoid radius R po iar increases to ~ 1.5, which is approximately the 
size of the asymptotic compactification half-period L/2 = tt/2. Thus it is clear from the prolateness 
and lack of deformation that the geometry around the symmetry axis is decompactifying. 

We now wish to characterise this decompactification. In figure El we plot a selection of black hole 
embeddings, now including the embedding of the symmetry axis to show its proper length. Note that 
we only show half of the full period, and thus reflecting the horizon and axis about X = generates 
the full compact period. The asymptotic compactification radius for half the period is L/2 = ir/2 here. 
We term the length of the axis for half the period, L ax i s , and it is given as, 

Laxis = / dpJe 13 ' 1 | x=7r / 2 (25) 
J po 

We see L ax i S decreases in these figures for larger black holes, but the horizon radii increase faster, 
resulting in an overall decompactification. In figure [7| we show the maximum value of the embedding 
coordinate X to contain half a period of both the horizon and the symmetry axis, 

Xmax = X{tt/2) + L axis (26) 

which is essentially the same as [Rpoiar ~t~ L ax i s ^j since the ellipse is such a good fit to the horizon 
geometry. We take this quantity to be the physically relevant (and coordinate invariant) measure of 
the compactification length near the axis. 
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Figure 4: Embedding of horizons into Euclidean space for a moderate black hole on the left, with 
po = —0.71, and the maximal black hole to the right with p = —0.18. Whilst we interpolate the 
numerical functions to determine the embedding, we plot here the actual positions of the lattice points 
for our highest resolution to give an indication of how the resolution varies with position on the horizon. 
For all the solutions, we may fit a prolate ellipsoid to the horizon, the dashed line, and in all cases find 
perfect agreement. 



These plots clearly show the axis decompactifying. Since we are only able to 'grow' black holes 
to po = —0.18, it is unclear whether; i) there exists a maximal mass black hole, or ii) whether the 
decompactification continues in such a way that arbitrarily large black holes may exist, eg. with 
Laxis —> or a constant in the limit of infinite mass. If there were a maximum mass, it seems likely 
from the figures that po can still be increased some way more before we would expect to reach it. 

3.2 Comparison with non-uniform strings 

We now compare these geometries with those of the critical uniform string (A = 0), and the most 
non-uniform strings (A ~ 4) found in ^2]> rescaling the asymptotic radius appropriately, and defining 
A as in that paper, 

X = U^-l) (27) 

where R ma x,min are the maximum and minimum horizon radii respectively. Whilst A = 3.9 was the 
most non-uniform string found there, the geometry and also thermodynamic quantities M, T, S appear 
to asymptote for large A (see also j2H]) and thus the A — > oo values are expected to be very similar 
to those at A ~ 4, probably only differing by a few percent. Strong evidence for this comes from the 
realization that a conical geometry forms at large A, as tested in [SO]. Once A is relatively large, say 
A ~ 2, only the geometry near the string 'waist' appears to change with increasing A as the cone forms. 
With the rescaling so that the asymptotic radius for one period of the solution is L = ir, we find the 
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Figure 5: On the left we plot the equatorial (red) and polar (blue) radii of the ellipses fitted to the 
horizon embeddings, against the parameter po specifying the size of the black hole. It is unclear what 
happens in the limit po — > where the p coordinate changes topology. It would be very interesting to 
know if R eq , R po i ar remain finite or not, and thus whether there is a maximum mass black hole or not. 
Since the geometry near the axis and horizon decompactifies, even though the largest black holes found 
have comparable radii to the corresponding asymptotic half period distance L/2 = tt/2, the horizons 
are still quite spherical. We plot the ratio of the radii in the right hand diagram. Again it is unclear 
what will occur in the limit po — > 0. 

following values, 



Critical string 



Rmax,X=0 
X m ax,X=0 

M x=0 
T x=0 



0.64 

7T 

2 

1.47 

0.250 

3.99 



(28) 



Highly non-uniform string R ma x,x=3.9 = 1-H 

X m ax, A=3.9 = 1>91 

M A=3 .9 = 3.38 
r A=3 . 9 = 0.184 
Sa=3.9 = 12-94 

Here R ma x,x is the maximal radius of the horizon, and X maX) \ is again the maximum value of X when 
embedding half a period of the strings into Euclidean space, using the metric (|22|) as for the black holes, 
and taking X = at the maximal radius of the horizon, so X = X max at the 'waist', the minimum 
radius. Now, there is no exposed symmetry axis, so X max is just the change in X when traversing the 
horizon for a half period. As noted in |2H] the proper distance (which is not equal to the embedding 
coordinate X) along the horizon increases with A indicating the geometry decompactifies there. 

The key observation of this paper is now evident from the previous plots of the embedded black 
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X 

Figure 6: We show both the embedding of the horizon, and also the symmetry axis for various 
solutions ranging up to the largest black hole found. These embeddings are for half a period, and 
should be reflected about X = to obtain the full geometry. The asymptotic radius of this half period 
is L/2 = it/2. We clearly see that the total space taken up by the solutions in the embedding X 
coordinate grows as the black hole increases in size, despite the symmetry axis decreasing in length, 
indicating the overall geometry near the horizon and axis decompactifies. We also see that even the 
largest black hole still appears to have considerable 'room' to increase its size still further. 




Figure 7: Plot showing X max , the 'length' taken up by half the horizon and symmetry axis when 
embedding into Euclidean space. This quantity offers a physical and coordinate invariant measure of 
the radius of compactification near the horizon. We see as the black hole grows, with increasing po, 
the axis decompactifies considerably. Again, it is unclear what happens in the po ~~ > limit. 
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hole geometry. Already for p ~ —0.35 we see the black hole equatorial radius R eq is equal to R m ax f° r 
the highly non-uniform A = 3.9 non-uniform string, which as stated above we take to be approximately 
equal to R ma x in the limit A = oo. And as po increases past this point, the black holes continue to 
become larger in radius. In addition we also see from X max that the geometry along the axis has 
decompactified as much as that of the very non-uniform strings for even quite small black holes with 
po — —1.0, and again the trend seems to continue for increasing po past this point. The implication 
is clear, that it seems difficult to imagine the black hole solutions making a transition at po = via a 
cone geometry to the most non-uniform strings, when they simply become 'bigger' than the most non- 
uniform strings already for po — —0.35. We may gain more insight into this result in figure[E]by plotting 
the embedding of the A = 3.9 non-uniform string, and the largest black hole relaxed (po — —0.18) into 
Euclidean space, including the symmetry axis of the black hole. Again we note that with increasing 
A past 3.9 we only expect the geometry in the cone region to change, and thus the geometry of the 
A = 3.9 string should be very close to that of the limiting string at A = oo (see |29] for curves of R m ax 
against A). 

Whilst we earlier claimed that resolution becomes limited near the horizon and axis for the large 
black holes, we find that the values of R eq , R po iar only vary by ~ 1% when doubling the resolution from 
70 * 210 to our highest resolution 140 * 420 for the solution with po — —0.35 that parallels the most 
non-uniform string horizon size. The length along the axis Laxis varies only a little more, around ^ 3%. 
Thus while decreasing axis resolution with increasing black hole size limits the ability of the algorithm 
to converge, the resolution is still high enough for the accuracy of our large black hole solutions to 
be high. For further comparison of quantities measured at different resolutions, see figure and the 
Appendix lAl 

3.3 Thermodynamics 

Now we turn to the thermodynamic quantities to see whether our comparison of the black hole/string 
horizon intrinsic geometry is parallelled in these independent observables. We might now reasonably 
expect the mass of the black holes to become greater than that of the most non-uniform strings, and 
the horizon temperature to become less. In figure El we plot the temperature, T, and mass M of the 
black hole solutions, now against the horizon entropy S. The mass is computed in two ways, firstly 
asymptotically from the metric [see equation (|18[) ]. and secondly by integration from the First Law. 
We clearly see very good agreement for these as expected. As discussed earlier in section f2.5| this is 
a good test of the elliptic equations, but does not test the constraints which are the really important 
quantities to check for this elliptic method, as they are not imposed directly. These constraints are 
tested explicitly in Appendix iBl 

On these plots we also show the same quantities for the non-uniform string branch up to A = 3.9. 
Again we emphasise that the A = 3.9 point probably lies very close to the A = oo point on this diagram. 
Our expectations are confirmed, with the temperature becoming lower, and the mass becoming higher 
than the most non-uniform strings already by po — —0.30, and the trend continuing for larger po, again 
reflecting the fact that the black holes become 'bigger' than the non-uniform strings. We also plot the 
behaviour of a 6-dimensional Schwarzschild solution, and find that since the axis is decompactifying, 
and as a consequence the black hole horizon geometry is only slowly deforming from being a sphere, 
the black holes closely reproduce this Schwarzschild behaviour. 

There is one further point we may observe from these plots. The non-uniform branch presumably 
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X 

Figure 8: An illustration of the intrinsic geometries of the maximum sized black hole found (po = —0.18) 
compared to a highly non-uniform string (A = 3.9) whose geometry is expected to be very close to the 
limiting A = oo solution, with a conical geometry near the waist of the horizon. Only half a period of 
the solutions is shown, and for the full geometry one should reflect about X = 0. Both are taken to 
have an asymptotic compactification half period L/2 = ir/2. We see immediately that the black hole 
is simply 'bigger' in all aspects. It has larger equatorial radius, and including the exposed symmetry 
axis, it takes up more 'room' in the X direction - ie. it has decompactified more. Indeed we also find 
it has a higher mass, and lower temperature. Furthermore, it looks very much as if we can increase the 
black hole size further still. The implication is that the A = oo string solution presumably cannot be 
connected (through a conical topology change) to this branch of black hole solutions. 
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terminates at A = do, very close to the A = 3.9 solutions plotted, so the whole non-uniform branch of 
solutions A = to oo then appears to lie very close to the black hole curves both in A4 and T against 
S. Since the branches do not appear to connect, we can think of no particular reason why this should 
be so, and presumably it is just an interesting coincidence. 

Finally we turn to the last (and truly higher dimensional) thermodynamic quantity, the binding 
energy n of our solutions. This turns out to be very small, and thus prone to numerical error. Referring 
back to (|21|) we recall n to be given by a ratio of asymptotic quantities, with numerator 3 62 — 3 C2 + 02- 
In figure E3 we plot the magnitudes of both a term in the numerator, and the numerator itself, and 
demonstrate that the numerator is relatively very small, ~ 10% of the terms making it up, indicating 
cancellations occur between the terms. This is a problem numerically since the asymptotic quantities 
a 2;^2) c 2 are already difficult to measure, and thus a quantities depending on detailed cancellations 
between them certainly should not be trusted. A further caveat is that for small black holes both the 
numerator and denominator are small, and their ratio is consequently extremely unreliable. Computing 
n from the asymptotics we find its value to be less than ~ 0.1 for these solutions, but the errors appear 
large, and we stress simply that it is small, and we do not feel we can give its value with certainty here. 
We may reassure ourselves that whilst n is 'noisy' in the numerical error, this is simply because n is 
close to vanishing, and Smarr's law is extremely well satisfied. Also in figure EH we plot the two sides 
of the Smarr relation (|20[). TS and (3 — n)M/4, against po. In addition, we also plot 3M/4, ie. setting 
n = 0, on the same plot, which lies so close to the curve with the actual measured n, that it is clear n 
is quantitatively small, and furthermore we cannot currently expect to measure it with any accuracy 
as argued above. Much higher precision would be required for this, and thus we leave determining the 
actual (small) values of n for future work. 

4 Discussion and Outlook 

We have shown that static black holes in pure 6-dimensional gravity compactified on a circle, ie. Kaluza- 
Klein black holes, may be found using elliptic numerical methods. As expected, for fixed asymptotic 
compactification radius, the small black holes behave as 6-dimensional Schwarzschild solutions. As 
they grow the geometry on the axis decompactifies relative to the fixed asymptotic radius, and the 
horizon deforms to a prolate ellipsoid. Since we are limited by numerical resolution, we are currently 
unable to probe whether the axis decompactifies indefinitely, and consequently black holes of any mass 
can be found, or whether instead there is an upper mass limit for the black holes. This is clearly an 
interesting issue to resolve in future work. 

The most interesting result we find is that whilst we are only able to compute black holes whose radii 
become approximately equal to the asymptotic compactification radius, these are already sufficiently 
large that they are simply 'bigger' than the most non-uniform string solutions constructed in ^2]; both 
in terms of horizon volume and mass. For the largest black hole solutions we found the horizons are 
still quite spherical since the axis geometry decompactifies, making 'room' for the horizon. It therefore 
appears that they should be able to further increase in size and mass, past the point where we are 
currently able to construct them, before any possible upper mass limit would be reached. Given these 
results, it therefore seems rather unlikely that this black hole branch of solutions can merge with the 
non-uniform strings via a conical geometry developing at the polar regions of the horizon, as suggested 
by Kol [SHj. This is despite the fact that we have excellent numerical evidence that the highly non- 
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Figure 9: Plot showing black hole temperature T and mass M against entropy S calculated by eqns. 
(|17|) and (|18[l. Three resolutions for the black holes are used, 36*106 in green, 70*210 in red, 140*420 in 
blue, and excellent consistency is found. In the figure for the mass, the dashed black line line gives the 
mass computed by integration from the First Law, showing excellent agreement with the asymptotically 
measured mass. The solid red lines give the behaviour of a 6-dimensional Schwarzschild solution, and 
for these black holes, despite their size becoming equivalent to the compactification radius, we see little 
deviation from this. We also show the same quantities plotted for the non-uniform strings as pink open 
squares. The dashed and solid straight lines show the limiting values for A = 0, and A — * oo. We 
include magnifications of the region where the non- uniform strings exist. Interestingly the non- uniform 
values lie very close to the black hole curves, although this appears simply to be coincidence. 
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Figure 10: Left; A plot of the magnitude of a term in the numerator of 'n', the binding energy, and the 
magnitude of the numerator itself. We see this numerator is ~ 10% of the value of one of its constituent 
terms, and thus detailed cancellations occur, and the value of 'n' is very small. Since the asymptotic 
quantities 02 , 62 , C2 are already difficult to compute accurately, we should not trust a numerical quantity 
that depends on detailed differences of these, and thus is very small in comparison. Right; A plot of the 
left- (TS as a solid black line) and right-hand ((3 — n)M/4) as blue points) sides of Smarr's formula, 
confirming that whilst n is very hard to determine accurately, Smarr's formula is very well satisfied by 
our data. To confirm that 'n' is indeed extremely small, we show the right-hand side computed with 
n = 0, ie. 3M/4 as red points, which gives a slightly worse fit, but only marginally. 

uniform strings do indeed exhibit the required conical geometry at their waist, which previously lent 
weight to this conjecture |50] . 

This is then very interesting geometrically, and raises the obvious question of whether the static 
non-uniform string solutions can be continued through the A — > 00 solution with its conical waist to 
a new branch of black hole solutions. Obviously while these would have horizons that do not wrap 
the extra coordinate, they would be distinct solutions from the ones we construct here, and at low 
mass (if they have a low mass limit) would presumably not look like the 6-dimensional Schwarzschild 
solution. 7 Similarly, if the black hole branch we partially construct here does turn out to have an 
upper mass limit, can this branch of solutions be continued through to a new string solution, distinct 
from that connected to the Gregory-Laflamme critical uniform string? We refrain from speculating on 
these questions (see the recent |28| for an interesting discussion of a variety of possibilities), deferring 
these issues until improved numerics can be performed that confirm the current results, and can extend 
the range of these elliptic methods so these questions may be tackled directly. We do make one further 
general comment here. As we have seen comparing the current work here with the previous work 
constructing the non-uniform string solutions ^5], the fact that the axis of symmetry is exposed for 
the black hole solutions completely changes the boundary conditions imposed on the problem. Thus, 
without proof, it would be dangerous to assume that continuation of a branch of solutions through 
topology change at a conical region must always be possible. 

7 For example the non-uniform string branch may join a branch of black holes that then undergoes additional topology 
changing, and so in the limit of small mass the horizons are Schwarzschild-like, but with more than one horizon per 
compactification period. 
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Whilst our computations were performed in 6-dimensions, in order to make contact with previous 
work constructing the non-uniform strings, it would be good to check the same behaviour occurs 
in 5-dimensions. Whilst the difference between 4-dimensions and more than 4 is very large, due to 
the additional curvature terms entering the Einstein equations from the rotation group of the axial 
symmetry, the difference between 5 and 6 dimensions is simply in coefficients entering these equations. 
Hence we would be very surprised if the 5 and 6-dimensional systems behaved qualitatively differently, 
but to be sure it would be good to check by constructing the non-uniform string and black hole solutions 
in 5-dimensions and comparing them. 

Our findings appear to be strongly related to the decompactification of the geometry near the hori- 
zon and axis. The reason the black holes become larger than the non-uniform strings is because the axis 
decompactifies making room for them. Therefore in order to make contact with realistic phenomenol- 
ogy, and thus really determine whether there are interesting strong gravity effects of compactification, 
it is clearly important to consider the problem again, but include some radius stabilising mechanism. 
Presumably once stabilisation is included, an upper mass limit for the black holes should be inevitable 
as the axis can't decompactify so easily. This is a sufficiently important phenomenological question 
that this should be checked explicitly, rather than just assumed, as if it turned out not to be true, or 
only be true for certain stabilisation mechanisms, this might provide new physical and observational 
constraints on compactifications that are totally independent of the familiar weak field constraints. It 
may also provide an important testing ground for the non-linear dynamics of stabilisation mechanisms. 
Additional matter, such as is required for stabilisation, simply adds elliptic equations to the problem, 
and no further constraints, and thus in principle can be easily incorporated. We note that, at least with 
weak stabilisation, the Gregory-Laflamme instability will occur as usual. Hence non-uniform string so- 
lutions will also exist, although of course it is then not obvious that they can be deformed to have a 
conical region of their horizon as in the unstabilised case. If they do behave in the same manner as 
for the unstabilised theory, and if, since the axis could not decompactify readily, the black holes are 
forced to have an upper mass limit at a lower mass than those probed in this paper, then possibly the 
topology change to the black hole branch that Kol suggests could occur after all. 
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A Appendix: Numerical details 

As in the previous works |18|I12|I19| we use a simple Gauss-Seidel method with second order differencing 
to relax the elliptic equations. The non-linear source terms are fixed, the resulting Poisson equations 
are relaxed, and then the sources are updated with the new solutions, and the boundary conditions are 
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refreshed. Repeating this, the elliptic equations are either completely relaxed and we find a solution, 
or convergence is lost at some point early in the relaxation, and all the metric functions diverge 
dramatically to nonsense. We impose the asymptotic boundary at finite p = p ma x which we typically 
take ~ 5, so several multiples of the half periodic compactification radius L/2 = n/2. In the next 
Appendix we show data for varying p ma x and demonstrate this has been taken large enough so as to 
be irrelevant. 

Essentially all our numerical problems come from the symmetry axis r = 0. Firstly, rather than 
discretising the grid in the p, x coordinates, generically one gains stability using p, £ with 

e = 4( X -vr/2) 2 , (29) 

since at x = 7r /2, all fields are even in (x — vr/2) and therefore linear in £. This was used successfully 
in the black hole on an RS brane to improve stability |19j . However, even with this modification the 
algorithm is horribly unstable, and even for the smallest black holes we find no convergence. The 
same problem was encountered in the earliest application of this method [18] . When dealing with 
a spherically symmetric scalar field <j) in polar coordinates, one is very familiar with terms such as 
^d r 4> in the equations of motion. Whilst for an elliptic relaxation these look as if they might destroy 
convergence, in reality they are not a severe problem, as long as the Neumann condition on ^ is imposed 
at r = 0. However, the ansatz (^Q) or Q generates more singular terms in the field equations. The 
exact form of our metric ansatz Q guarantees that only one of the elliptic equations is effected, but 
we find that in the equation for C, 

jAC=^(e- 10C / 3 -l)+... (30) 

where J A = J (d 2 + = (d 2 + d 2 .). The remaining terms have the more usual 1/r multiplying 
derivatives. Obviously the above term is finite as C ~ r 2 near the axis. However since we do not 
impose that C goes quadratically near the axis, and instead it emerges from a combination of the 
elliptic equations and the constraints, during the early stages of the relaxation this term generically 
destroys convergence. 

We deal with this term as in JH]- The second derivative terms in the constraint equation G p x 
are simply C px and thus have characteristics compatible with integrating C over the (p, x) domain. 
However this equation has no such singular term as that above, and any solution for C integrated away 
from the r = axis has very good quadratic behaviour in r near there. 8 Thus using this constraint, we 
integrate for C, but call this function C2. Since it has very good properties near the r axis, we calculate 
the one singular term in the elliptic equation for C using this function C2, rather than C. Whilst this 
seems circular in nature, and we offer no proof why this should converge so that C = C2 finally, in 
practice this does indeed happen, and the method becomes very stable. Given the characteristics we 
need two initial data surfaces, one at constant x, the other at constant p. For x we fix C2 = C at 
X = tt/2, which includes the r = axis for p < 0, and this ensures the quality of the C2 behaviour is 
good near the symmetry axis. For the other initial surface we take the horizon, and fix C2 using the 
condition (|13|) (but now for C2) that the horizon temperature is a constant, giving 

C2( P0 ,x) = (B-A)-(B- A)\ x=n/2 + I log (jj^) ■ (31) 

8 One might imagine using just G p x and 2 elliptic equations for A,B, but we have been unable to find a scheme like 
this that worked in practice due to the non-local nature of the integration for C. 
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Figure 11: Plots of the metric functions C and C2 for po = —0.71. Main error of C2 originates from 
the coordinate singularity at p = 0, x = 7r /2, but its effect for C through the elliptic equation is small, 
and the difference of their magnitude is very small compared to C. 

As discussed earlier, this condition is not used when solving the elliptic equations, where we instead 
use the (G p p — G x ) constraint. However, it gives a more stable initial boundary condition for the C2 
integration than simply setting C2 = C there directly (which destroys convergence). We then integrate 
C2 from these two boundaries by quadrature, 

C2 = f P dp [ X d X F(p, X ) + C(p , X) + C(p, vr/2) , (32) 

J po J tt/2 

where F is the 'source' term in the constraint equation. 

The reason this method to compute the singular term in the elliptic equation for C appears to 
work is that in practice the contribution ( e - 10C V 3 — l)/r 2 is only significant near the r = axis, and 
away from there this source term dies away more quickly than the other terms as it is suppressed by 
1/r 2 . Thus whilst the process appears very non-local, involving integration over the lattice during 
the relaxation, which is not very 'gentle' and might destroy convergence, actually it only has an effect 
localised at the axis. Furthermore the function C is generically much smaller than the other metric 
functions A, B, so C appears to have relatively little effect on the solution anyway. 

Since we do find stable converged solutions we may check that C is equal to C2, and indeed 
comparing these globally over our domain is an excellent check that the constraint equations are 
enforced. In figure^2 we show C and C2 for a black hole with p^ = —0.71, and we see the difference of 
the functions is very small compared to C, and furthermore C is very small compared to A,B. Thus 
this method, whilst appearing rather mysterious, does seem to give very good results in practice. 

We reconstruct the values of A, B at the coordinate singularity p = 0, x = 7r /2 (C of course is zero), 
and since the p, x coordinate system is singular there, it is easiest to consider this in the non-singular 
r, z coordinates. Then from the axial symmetry at r = and the reflection symmetry at z = tt/2, the 
metric functions have expansions, 

A(r, z) = k + hr 2 + k 2 (tt/2 - zf + 0(r 4 , (vr/2 - z)\ r 2 (vr/2 - z) 2 ) (33) 
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Figure 12: Plots of the weighted constraint equations for pq = —0.71. 

and similarly for B, and thus using the relation of r, z to p,x, we may simply compute how to interpolate 
the values at the coordinate singularity r = 0, z = tt/2 from the neighbouring points in the p, \ grid. 



B Appendix: Numerical checks 

We find second order scaling in all physical quantities such as the mass, temperature and entropy. 
However, the resolutions are sufficiently high that increases yield very little change in the quantities. 
Observe the earlier figure where we have plotted the thermodynamic quantities for 3 different resolu- 
tions, which give extremely similar results where multiple resolutions may be relaxed. This gives much 
confidence that the accuracy of these solutions is high. Certainly our main conclusion is seen in this 
figure, that the black holes become larger in mass and entropy than the most non-uniform strings, and 
we see this is totally unaffected by changes in resolution. 

As discussed in the main text, it is essential to explicitly test that the constraint equations are well 
satisfied, as these are not imposed directly, but only via boundary conditions and the CR relations. For 
po = —0.71 we plot in figure^Jthe weighted constraints \P and <I>. We see that they are suitably small, 
rising to their maximum near the symmetry axis or coordinate singularity. However, it is very difficult 
to interpret these constraint violation values in terms of their physical effect. A nice check that these 
small violations are sufficiently small that the physics of the solution is unaffected by them is given in 
Table Q where we show the average values of the weighted constraints and C2 — C (which also gives 
a measure of how well G p x is satisfied) over p and x f° r three resolutions with different pq. As the 
numerical resolution is increased, the averaged constraint violation values decrease significantly (not 
quite as quickly as second order scaling, but then our discretisation geometry is rather complicated so 
this would not be expected), indicating the constraints become increasingly well satisfied, as we would 
hope for. The geometry and other properties of the solutions varies very little as the resolution is 
increased, and thus the constraint violations must be very small in terms of their physical effect. 

In order to assess the absolute physical error in these small constraint violations globally, we advo- 
cate comparing the values of C and C2 shown in the previous Appendix over the whole domain. Since 
these agree extremely well, this is again excellent evidence that G p x is effectively very well satisfied, 
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Table 1: This table shows averaged violations of the weighted constraint equations and C2 — C for 
three resolutions and three different black holes. The average of absolute values is taken over the whole 
domain. 
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\SS/S\ 


\6M/M\ 


Pmax — 4.6 


3.5 x 10" 4 


1.3 x 10~ 3 


1.1 x 10~ 2 


2.5 


3.6 x 10" 3 


1.9 x 10" 2 


9.3 x 10" 2 


1.4 


5.1 x 10~ 2 


1.8 x 10" 1 


4.0 x 10" 1 



Table 2: This table shows the variations of the temperature, entropy and asymptotic mass for a black 
hole with po = —0.36 using different p m ax- The variation of the quantities is defined as 5T/T = 
1 - T/T Pmax=5 . 7 , dividing by T for p max = 5.7. 

since C2 is integrated from this constraint, and C is obviously derived from the elliptic equations. 
Yet another physical test we may perform is to compute the horizon temperature as a function of x 
along the horizon. Again the G p x constraint should ensure this is constant, yet it is the (G p p — G\) 
constraint that we actually impose at the horizon for the elliptic equations. Thus if G p x is well satis- 
fied, then the CR constraint structure is working well. In figure El we plot the maximum variation, 
5T error = (T max — Tmin) /T, of the temperature for varying po for our highest resolution, and two lower 
resolutions. These variations are small, implying the constraints are indeed well satisfied. Furthermore 
they decrease very nicely with increasing resolution indicating the constraints are behaving well nu- 
merically, and are free from systematic violations. We find maximum variations for the largest black 
holes, as expected as gradients build up near the horizon at x = 7r/2 due to the limited resolution at 
the axis. However the variation is still only ~ 1% for the largest black hole we relaxed. 

In the TableElwe show the temperature, entropy and mass for a black hole using different p ma x- We 
see these quantities (and indeed all others) hardly change, indicating our choice of p ma x is sufficiently 
large. 
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Figure 13: Plot showing the variation of the temperature on the horizon, 5T error = (T max — T min ) /T, 
for three resolutions; our maximum 140*420 and two lower resolutions, 70*210 and 36*106. 

C Appendix: Demonstration that First Law does not test the con- 
straints 

In this brief Appendix we demonstrate our claim that deriving the asymptotic mass by integrating 
the First Law dM = TdS along our branch of numerical solutions (for fixed asymptotic radius) does 
not test whether the constraint equations are satisfied. This is also true for the related Smarr law. 
Thus while it is useful to check that the First Law is satisfied, as it checks the elliptic equations are 
well satisfied and have boundary conditions imposed compatibly with their regular singular behaviour, 
we should not be lulled into a false sense of security. We satisfy the elliptic equations directly in this 
relaxation method, and to high accuracy, whereas the constraint equations are imposed indirectly, via 
the boundary conditions. Hence these are the equations we should worry may have numerical errors, 
and it is crucial to separately check these, as was done in the previous Appendix. 

The First Law can be classically derived, eg. as in j^E]- Here we simply sketch the derivation, con- 
sidering which components of the Einstein tensor are involved, and therefore which can be numerically 
tested by the First Law. Consider the expression, 



for a manifold Ai with metric g^ v . This is not the action, since for the action we must subtract the 
Gibbons-Hawking term at the boundaries 5M. Clearly S (unlike the true action) vanishes for any 
solution of the equations of motion, as the Ricci scalar will always vanish locally. 

The first law can be derived in our static case by considering g^ u = gffi + Sg^u, where gj£) is a 
static solution of the Einstein equations, and the perturbation Sg^ u also satisfies the static linearised 
perturbation equations. Thus, S vanishes when evaluated on both g^ and g. However, in the usual 
way we can write, 



where d n is the derivative normal to the boundary, and V, W give boundary terms that arise to eliminate 
derivatives of Sg in the integral term. Since the first 3 terms all vanish, we are left simply with the 
boundary terms, and these give rise to the First Law, linearised in 5g, when evaluated on the horizon 




(34) 



S(jg a p) ~ S(g { V) = [ ffxyJ^G^ig^SsT + W^gT + W^dJg^sM + °( 6 9 2 ) (35) 



JM 
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and asymptotically. However, let us now consider the above in our numerical context. If numerically we 
see the First is well satisfied, does this imply all the Einstein equations, both elliptic and constraints, 
are therefore well satisfied? 

Naively this appears so. From the terms ^(g^), S(gffi) we test the weighted average of the Ricci 
scalars, R(g) and R(g^), and from the linear variation term we test a weighted average of G^i^g^). 
However, given that our metric (^Q) is diagonal, and further more, our perturbation 5g is therefore also 
diagonal, these three terms then only test the diagonal components of the Einstein tensor. Thus these 
weighted integrals simply do not involve the G r z Einstein equation (or equivalently G p x ). Even worse, 
due to the conformal invariance of the r, z (or equivalently /, g) block of the metric, the perturbation 
is restricted so 5g rr = Sg zz and therefore the integrals also do not involve (G r r — G z z ) (or equivalently 
(G p P -G\)). 

Now it is clear that if we take any solution of the elliptic equations, g(°\ totally ignoring the 
constraints, and perturb this by 5g which again only satisfies the static linearised elliptic equations, we 
can perform exactly the same manipulations to obtain the above equation, and consequently the usual 
First Law. Hence for 'solutions' only obeying the elliptic equations we would still see the First Law 
being well observed, even though we had made no attempt to satisfy the constraint equations. 

Thus, our specific form of the metric exactly ensures the Einstein equations that are the 'constraints' 
in this elliptic context, G r z and (G r r — G z z ), do not affect the First Law, even if they are not satisfied due 
to numerical error. For exactly the same reasons, the closely related Smarr law also has no dependence 
on the constraint equations, when using our metric choice, and similarly cannot provide a numerical 
test of the constraints. 
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